Quasiparticle scattering and local density of states in graphite 
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We determine the effect of quasiparticle interference on the spatial variations of the local density 
of states (LDOS) in graphite in the neighborhood of an isolated impurity. A number of characteristic 
behaviors of interference are identified in the Fourier transformed spectrum. A comparison between 
our results and scanning tunneling microscopy (STM) experiments could provide a critical test 
of the range (of energy) of applicability of the Fermi liquid description of graphite, where some 
evidence of the breakdown of Fermi liquid theory has recently been discussed. Moreover, given 
the similarity between the band structures of graphite and that of nodal quasiparticles in a d-wave 
superconductor, a comparison between results in the two materials is useful for understanding the 
physics of the cuprates. 
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PACS numbers: 68.37.Ef, 81.05.Uw, 71.10.Ca 

There has been considerable interest recently in 
Fourier transform scanning tunneling spectroscopy (FT- 
STS) measurements, which may reveal some rather in- 
teresting electronic properties of conducting materials. 
In particular, recent experiments performed on certain 
cuprate high Tc superconductor a^'^'^'^'^'^'^ , promise to 
shine some light on the low energy electronic excitations 
of the superconducting state, and on the pseudogap. 

In this paper we study a simple two dimensional sys- 
tem, graphite, which is generally believed to be a Fermi 
liquid. (Note, however, that this belief has been ques- 
tioned in Refs. 0.) Graphite has many similarities to 
the cuprates; its band structure consists of two bands 
that touch only at the six corners of the Brillouin zone 
(BZ). As a result, graphite has nodal linear dispersing 
quasiparticle excitations - analogous to those of the Dirac 
equation. This makes it in some ways analogous to a 
d-wave superconductor, where the nodal quasi particles 
are expected to have a (highly anisotropic) Dirac spec- 
trum. Thus, our interest is in obtaining an experimen- 
tally testable understanding of the implications of quasi- 
particle interference, both as a test of the range of validity 
of Fermi liquid theory in graphite, and to provide a fidu- 
ciary point for discussions of the more complex situation 
in the cuprates. 

We use a T-matrix approach^'^"'^^ to obtain the quasi- 
particle interference spectra for various energies, in the 
presence of a single impurity. Our predictions can be 
easily tested experimentally, as very clean graphite sam- 
ples are relatively easy to obtain. The resulting FT-STS 
maps can be quite complex, and contain regions of high 
intensity. Depending on the energy, these regions can 
be circular, triangular, or hexagonal, as shown in Fig. 0] 
For example, at low energy, the dominant features in 
the FT-STS maps are circles centered about the corners 
and the center of the BZ. The contours evolve with in- 
creasing energy: their radii increase, the circles centered 
about the corners become triangular, while the circle in 
the center becomes a hexagon. Beyond the critical energy 
defined by the Van Hove singularity, the topology of the 



contours changes, and the LDOS exhibits hexagonally 
shaped lines of high intensity centered about the cen- 
ter of the BZ. With increasing energy these lines become 
smaller and circular again, and they disappear altogether 
at even higher energies. 

We now review some general properties of graphite and 
explain our T-matrix calculation of the FT-STS maps in 
graphite in the presence of a single impurity. 

Graphite has a hexagonal lattice with one valence elec- 
tron per atom and two electrons per unit cell. The band 
structure of graphite consists of two bands with an en- 
ergy dispersion approximatediSi by the simple particle- 
hole symmetric tight-binding band structure: 



e±(fc) 



(1) 
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where a is the distance between nearest neighbors, and 
7 Ri leV. For simplicity we will set 7 = leV and a = 1 
in the rest of the calculation. The two bands touch at 
the corners of the BZ, at points Kq = 27r/9(±2A/3, 0), 
27r/9(±\/3, ±3). The band structure for graphite is de- 
picted in Fig.^ At low energies, for wavevectors close to 
the corners of the BZ, the dispersion of the nodal quasi- 
particle excitations is linear e±{k) ~ ±37|fc — Ko\/2. 

The single particle density of states (DOS) of graphite 
can be easily obtained by integrating the spectral func- 
tion, po{uj) = / (Pk/{4TT'^)A{k,uj), where the integral is 
performed over the BZ, and 



A(fc,w) = -2 ^ Im|[tj-ej(fc)-|-i(5] ^| 



(2) 



In the numerical evaluation of the single particle den- 
sity of states we took a finite quasiparticle inverse life- 
time 6 = 0.02eV. The result is sketched in Fig. |21 The 
tiny DOS at = 0, as well as the rounding of the 
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FIG. 1: The equal energy contours \ei{kx, ky)\ — \e2{kx, ky)\ 
in graphite. The six points marked by dots correspond to the 
corners of the BZ (also the points with ui — 0), and the BZ 
is indicated by the dashed lines. The band structure at small 
energies consists of circular energy contours centered about 
the corners and the center of the BZ. At larger energies, the 
circles change shape and size, and above some energy u — leV 
they become hexagons and then circles centered about the 
center of the BZ, to disappear altogether for energies larger 
than Lo = 3eV. 



van Hove singularities, are due to the finite quasipar- 
ticle lifetime. When the quasiparticle inverse lifetime S 
approaches zero, the density of states at the Van Hove 
singularity diverges Pq{wvh) oc ln(l/(5), and the density 
of states at zero energy goes to zero, po{0) oc S\n{l/S). 
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FIG. 2: The single particle density of states p{uj) plotted as a 
function of energy lj. We note the v-shaped DOS of small en- 
ergy and the presence of the van-Hove singularities at energies 
of ±leV. At energies \uj\ > 3eV the DOS goes to zero. 

We note that the density of states is symmetric about 
the origin ui = and is v-shaped for energies < leV. 
At the points ujvh = ±leV, the density of states shows 
van Hove singularities. This is the point where equal en- 
ergy contours change topology in the band structure (see 



Fig. P). As expected, the DOS goes to zero at energies 
\uj\ > ±3eV. 

Following Refs. |3 and|2 we compute the effect of a 
single impurity on the LDOS in graphite. We define a 
finite temperature (imaginary time) Green's function, 

G(fci , A:2 , r) = -Tr e-'^^^-^^) T. (r) (0) , (3) 

where K = H - fiN, e"'^" = Tre"''^, and T^ is the 
imaginary time ordering operator. The impurity scat- 
tering problem can be solved by computing the Fourier 
transform of the Green's function from the T-matrix 
formulationiSiii: 

G{ki,k2,iuJn) = Go{ki,iUn)T{ki,k2,iuj„)Go{k2,iuj„),{4) 
where 



Goik,iuJn) = 
and 











(5) 



T{ki,k2,iuJn) = V{ki,k2) 

+ Viki,k')Go{k', l^n)T{k', k2,liUn). (6) 

k' 

We will assume that the impurity scattering potential 
is very close to a delta function so that V is independent 
of k and k' , 



V = v 



1 1 
1 1 



For this case we can solve Eq. ©, and obtain 

T(zc^„) = [I-Vd^Goik,iu;n)]-'V 

V 



(7) 



(8) 



where / is the 2x2 identity matrix, a — 3^/3/2 is the 
area of the real space unit cell of the honeycomb lattice, 
and the integral over k is performed over the entire Bril- 
louin zone. In the neighborhood of the impurity, spatial 
oscillations of the local density of states are induced. It 
is easier to extract information about the quasiparticle 
characteristics from the Fourier transform of the LDOS: 



p{q,uj) ^iYg{k,q,uj), 



(9) 



where g{k,q,u;) = J2^=l,2(^idk,k+q,UJ)-G*,{k+q,k,uJ), 
and G(fc, k + q,uj) is obtained by analytical continuation 
iujn Lo + iS from G{k, k + q, iun)- We obtain 



Piq,u.)^-2j — 
4=1,2 



(10) 



r(^) 



Lu — ei{k) + i6][Lu ~ ei{k + q) + i6] 
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where T{uj) = v/{l - vd J d^k/{ATr^)[Gmk,uj + iS) + 
GQ^(fc, w + ii5)]}. Ea. (|ll|l is analyzed numerically and the 
resulting LDOS is plotted in Fig. 01 For our numerical 
analysis we pick values of w = 2eV and 6 — 0.02eV. 
Also we perform our numerical integration on a 288 x 240 
grid. We plot our results on a 144 x 120 grid. We have 
confirmed that the results have converged in the sense 
that changing the size of the grid leads to no observable 
changes in the results. 

At low energies, our analysis reveals that the dominant 
features in the LDOS spectra are circular contours of ra- 
dius iuj/Sj, centered about the corners and the center of 
the BZ. The shape of the contours changes with energy, 
and their size increases. At the energy uj = ujvh = leV, 
the spectral features change to hexagonal shaped con- 
tours centered about the center of the BZ. With increas- 
ing energy the contours become smaller and circular, and 
disappear for energies above w = 3eV. 

A few comments about our results are in order. We 
note that the intensity inside and outside the circles is not 
zero but finiteii, which, depending on the quasiparticle 
inverse lifetime (5, can blur the interiors of the circles, 
especially at low energies. 

The topology of the contours is identical for negative 
and positive energies. However the weight, and in some 
cases the sign of the peaks relative to the background, 
as well as the average spectral intensities are different, 
in particular at low energies. The asymmetry of the FT- 
STS spectra with respect to energy is not a direct band 
structure effect, since for graphite the band structure is 
particle-hole symmetric. It is an effect of the impurity 
scattering. We will discuss this asymmetry also in con- 
nection to the effect of the impurity scattering on the 
average density of states in Fig. O 

For positive energies, the results presented here do not 
depend sensitively on the strength of the impurity poten- 
tial, V. We have repeated our calculation for the FT-STS 
spectra for different values of v, including the limit of 
infinite scattering strength ("unitarity limit") and the 
changes are not qualitatively significant. However, at 
negative energies, our FT-STS pictures evolve with in- 
creasing V, and at large scattering strengths they become 
similar to their positive energy counterparts. 

Fermi liquid theory is most robust at small probing en- 
ergies, and at these energies we would expect the exper- 
imentally observed features to have the best agreement 
with the theory. As noted in Fig. ^ at small positive en- 
ergies the amplitude of the LDOS oscillations due to the 
impurity scattering is lower than at high energies, which 
may make the experimental observation of these oscilla- 
tions more difficult. However, the fractional amplitude of 
the LDOS oscillations at small negative energies is larger 
than at higher energies. 

The Van-Hove singularities which appear at high ener- 
gies, are a "big" feature, so survive in some form unless 
there is a very serious failure of the quasiparticle pic- 
ture. However, it is still possible that a strongly energy 
dependent quasiparticle lifetime could considerably mute 



the Van Hove singularities, while leaving the Fermi liquid 
behavior intact at much lower energies. 

We note that at low energy the position of the 
dots/circles is determined by symmetry, while at higher 
energies the shape and the position of the high intensity 
lines is sensitive to the assumed band structure. Another 
interesting point is that in a two dimensional system, 
the general structure of singularities is lines instead of 
points^"^. This yields lines of high intensity regions in the 
in FT-STS spectra due to quasiparticle interference, as 
opposed to pointsi^. 
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FIG. 3: The spatially averaged single particle density of states 
p(cj), in the presence of an impurity of strength a) i; = 2.5eV, 
and h) V — oo. 

Another interesting aspect we analyzed is how the spa- 
tially averaged single particle density of states p{uj) = 
Pa{uj) + 2cimpp{q = 0,0;) changes with the scattering 
strength v. Here Cimp is the concentration of impurities. 
We assume that Cimp is small, so the effects of the im- 
purities are uncorrelated. The resulting density of states 
depends on the value of the impurity potential. For com- 
parison, in Fig. 13 we give p{uj) for v = 2.5eV (a), and for 
V = oo (b). We fix the impurity concentration to c = 1%. 
As expected, in the unitarity limit the effect of the im- 
purity is to create a sharp peak centered at zero energy. 
For smaller impurity potential, the effect of the impurity 
is to create a smaller and broader peak at negative en- 
ergy. One also notes the existence of a sharp anti-bound 
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state at an energy larger than 3eV. The position of the 
anti-bound state moves to larger and larger energies as 
one increases the scattering strength. 

We also performed an analytical calculation of the lo- 
cal density of states as a function of position at low en- 
ergy, where we can approximate the spectrum by a linear 
(Dirac) dispersion lo = ±tii?|A:|. We focused on the case 
of a single impurity scattering. For points far from the 
impurity we found the correction to the density of states 
due to impurity scattering to have the form: 




where vf = 37/2 is the Fermi velocity. Here m = 2 
is the number of the independent inter-nodal scattering 
vectors (that cannot be obtained from other scattering 
vectors by translation by a reciprocal lattice vector or 
by reflection). We can take, for example, Qo = (0,0), 
Qi = 27r/9(2V3,0), Q2 = 27r/9(%/3, 3). The coefficients 
Ai and 4>i depend on energy and on the impurity charac- 
teristics. As expected, at large |a;|, the density of states 



decays as If one takes into account also the finite 

quasiparticle lifetime the above equation for the density 
of states changes to 

p{x,u:)^ p{x,u:)e-\^\^l--, (12) 

where r = t{ijj) is the quasiparticle lifetime. A measure- 
ment of the density of states as a function of position 
thus reflects the position of the nodes {Qi), the quasipar- 
ticle dispersion, vp, and the energy dependence of the 
quasiparticle lifetime. 
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FIG. 4: The FT-STS quasiparticle interference maps in graphite in the presence of a single impurity. We plot —p{q,Lo) at the 
energies uj = 0.17, 0.4, 0.8, 1, 1.1, 2, —0.17, —0.8 and — l.leV . The x and y axis of each plot are qx and qy, and we display our 
results for — 1.157r < ga, < I.IStt and — tt < Qy < "k. For clarity we draw the BZ in dashed lines. Note the change in the contours 
topology with increasing energy. The plots are drawn using a normalized linear gray scale, with white being the highest (1) 
intensity, and black being the lowest (0) intensity, as indicated. The actual values of the FT-LDOS corresponding to the lowest 
(0) and highest (1) intensity arc different for each energy, and we find them to be (—0.12, 1.73) for 0.17eV, (—0.04, 2.16) for 
0.4eV, (0.27, 3.6) for O.SeV, (0.45, 14.3) for l.OeV, (0.41, 3.96) for l.leV, (0.22, 6.45) for 2.0eV, (-8.2, -3.26) for -0.17eV, 
(-1.1, 4.07) for -O.SeV, and (0.29, 5.22) for -2.0eV in arbitrary units 



